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In  two  previous  Air  Force  Office  of  Scientific  Research  (AFOSR)  studies. 
References  1  and  2,  McDonnell  Douglas  Corporation  investigated  the  dynamic 
response  of  an  aerosurface.  Figure  1.  in  a  subsonic  airstreom.  Structural 
nonlinearities,  in  the  form  of  discrete  nonlinear  root  support  springs  Kq 
and/or  were  present  in  the  system.  These  studies  addressed  two  types 
of  nonlinear  springs.  Figure  2,  which  are  characteristic  of  those  encountered 
in  missile  systems  with  freeplay  or  deadband  in  the  aerosurface  root  support. 

The  present  research  program  is  expanding  on  these  previous  studies  for  an 
aerosurface  with  discrete  structural  nonlinearities.  The  nonlinear  system 
dynamic  behavior  is  being  investigated  in  detail  through  analytical  and 
experimental  activities.  The  research  activity  is  led  by  McDonnell  Douglas 
Missile  Systems  Company  (MDMSC)  with  support  from  the  University  of 
Missouri-Rolla  (UMR)  and  overall  program  direction  from  AFOSR.  Nonlinear 
aeroelastic  system  dynamic  analyses  are  being  conducted  by  MDMSC.  Design, 
fabrication,  and  testing  of  aerosurface  models  containing  discrete  structural 
nonlinearities  is  being  done  at  UMR. 

During  the  previous  work,  as  discussed  in  References  2  and  3,  it  became 
apparent  that  there  existed  regions  of  flow  dynamic  pressure  where  the  nature 
of  the  aerosurface  dynamic  response  was  not  understood.  These  regions  were 
most  apparent  for  an  aerosurface  containing  root  spring  freeplay  plus  preload 
nonlinearities.  Figure  2(b).  These  regions  appeared  in  a  limited  range  of 
dynamic  pressure  as  illustrated  by  the  data  presented  in  Figure  3. 
Additionally,  freguency  analyses.  Reference  2,  of  the  simulation  results  in 
these  regions  indicate  that  the  nonlinear  system  time  history  response  is 
comprised  of  several  frequencies  which  are  not  multiples  of  each  other. 

Analyses  performed  with  McDonnell  Douglas  funding,  since  the  previous 
AFOSR  studies,  indicated  that  the  response  in  these  regions  is  probably 
chaotic  in  nature.  Phase  plane  plots  of  root  pitch  response  for  dynamic 
pressures  within  and  on  either  side  of  the  region  of  unexpected  dynamic 
response  are  shown  in  Figure  4.  The  changing  nature  of  these  phase  plane 
plots  over  this  range  of  dynamic  pressures  is  characteristic  of  chaotic 
motion.  Reference  4.  Potential  existence  of  chaos  in  a  relatively  simple 
dynamic  system,  such  as  the  aerosurface  configuration  studied  in  these 
previous  research  programs,  is  very  significant. 

The  present  research  program  is  expanding  on  these  previous  results.  The 
dynamics  of  the  nonlinear  system  shown  in  Figure  5  is  being  investigated  both 
analytically  and  experimentally.  This  configuration  is  slightly  different 
from  that  shown  in  Figure  1.  This  system  has  root  support  pitch  (h)  and 
plunge  (a);  degrees  of  freedom  and  is  consistent  with  the  classical  two 
degree  of  freedom  flutter  model  such  as  discussed  in  Reference  5. 

During  the  completed  Basic  Year,  analytical  tools  have  been  developed  to 
investigate  the  time  domain  response  of  a  rigid  aerosurface  having  discrete 
nonlinearities  in  both  root  support  degrees  of  freedom.  These  analytical 
methods  are  applicable  to  either  steady  or  unsteady  aerodynamic  forcing 
function  representations.  An  aerosurface  root  support  mechanism  has  been 
designed  and  fabricated.  This  mechanism  provides  discrete,  independent  and 
adjustable  nonlinear  stiffness  for  the  two  aerosurface  root  support  degrees  of 
5470D/4 
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freedom.  A  rigid  aerosurface  has  been  designed  and  fabricated.  In  this 
context,  "rigid"  indicates  that  the  natural  frequencies  of  the  cantilevered 
aerosurface  are  well  above  the  modes  of  the  combined  aerosurface  and  root 
support  springs.  Wind  tunnel  testing  of  the  rigid  aerosurface  and  nonlinear 
support  mechanism  combination  has  been  initiated.  Testing  has  identified 
linear  system  flutter  points  and  demonstrated  nonlinear  system  limit  cycle 
responses  for  a  number  of  root  support  stiffness  combinations. 

During  the  upcoming  second  program  year.  Option  I,  wind  tunnel  testing 
will  be  completed  for  the  rigid  aerosurface.  Design  of  a  flexible  aerosurface 
will  be  completed  and  fabrication  begun.  Detailed  analytical  studies,  for 
both  the  rigid  and  flexible  aerosurface,  will  be  performed  during  the  coming 
year.  In  addition,  dynamic  testing  of  the  rigid  aerosurface  is  to  be 
initiated.  Test  to  analysis  correlation  studies  will  be  completed  for  data 
obtained  with  the  rigid  aerosurface.  These  analytical  and  experimental 
studies  will  be  directed  toward  identifying  potential  regions  of  chaotic 
response  in  the  nonlinear  aeroelastic  system. 

A  description  of  the  overall  research  program  is  presented  in  the 
following  section.  Detailed  discussions  of  the  Basic  Year  accomplishments  are 
provided  in  Section  3.  This  is  followed  by  discussions  of  planned  research 
for  the  second  or  Option  I  year.  Section  4;  conclusions  reached  to  this  point 
in  the  program.  Section  5;  a  listing  of  research  personnel  and  interactions. 
Section  6;  and  applicable  references.  Section  7. 


2 . 0  PROGRAM  SUMMARY 

This  is  a  three  year  research  program  to  investigate  chaotic  response 
behavior  Cor  aerosurfaces  containing  discrete  structural  nonlinearities.  The 
program  is  an  extension  of  previous  AFOSR  studies  and  contains  both  analytical 
and  experimental  activities. 

The  dynamic  behavior  of  a  rigid  aerosurface  has  been  investigated 
analytically  and  experimentally  during  the  first  year,  or  Basic  Year.  The 
rigid  surface  test  and  analysis  activities  are  to  be  continued  during  the 
second,  or  Option  I  year.  In  addition,  a  flexible  aerosurface  will  be 
designed  and  fabrication  begun  during  the  second  year.  The  third  or  Option  II 
year  of  the  program  will  move  to  test  and  analysis  of  the  flexible 
aerosurface.  A  description  of  the  planned  activities  for  each  of  these  years 
follows.  This  is  followed  by  a  description  of  the  approach  being  taken  to 
accomplish  these  planned  activities. 

2.1  Program  Plan 

The  activities  performed,  or  to  be  performed,  during  each  year  of  the 
research  program  are  described  in  the  following  paragraphs.  The  overall 
objective  through  out  the  program  is  to  obtain  an  increased  understanding  of 
the  dynamic  response  phenomenon  for  aerosurfaces  with  discrete  structural 
nonlinearities.  The  research  is  extending  the  fundamental  understanding  of 
nonlinear  system  dynamics  and  the  potential  characterization  of  this  motion 
within  the  context  of  chaotic  behavior. 

Basic  Year  -  Analytical  studies  have  been  performed  for  a  range  of  rigid 
aerosurface  configurations  and  various  root  spring  stiffnesses  and 
nonlinearities.  Test  apparatus  has  been  designed  and  fabricated  to 
experimentally  demonstrate  the  nonlinear  behavior  of  a  rigid  aerosurface 
containing  discrete  structural  nonlinearities.  Wind  tunnel  testing  for  the 
rigid  aerosurface  conf iguration  has  been  initiated  and  evaluation  of  the 
results  of  the  wind  tunnel  tests  is  underway.  Initial  design  and  fabrication 
of  the  rigid  aerosurface  dynamic  test  setup  was  completed  during  the  Basic 
Year . 

Option  I  Year  -  The  wind  tunnel  testing  for  the  rigid  aerosurface 
configuration  will  be  completed.  A  mechanically  excited  "feedback"  system 
will  be  used  to  perform  further  dynamic  testing  of  the  rigid  aerosurface 
configuration.  The  results  of  these  wind  tunnel  and  dynamic  tests  will  be 
evaluated.  Included  in  this  evaluation  will  be  test  to  analysis  correlation 
studies.  Analytical  studies  for  a  range  of  flexible  aerosurface 
configurations  and  various  root  spring  stiffnesses  and  nonlinearities  will  be 
performed.  The  design  will  be  completed  and  fabrication  initiated  for  test 
apparatus  to  experimental ly  demonstrate  the  nonlinear  behavior  of  a  flexible 
aerosurface  containing  discrete  structural  nonlinearities. 

Option  II  Year  -  A  mechanically  excited  "feedback"  system  will  be  used  to  test 
the  flexible  aerosurface  configuration.  Wind  tunnel  testing  will  be  performed 
for  the  flexible  aerosurfaces  configuration  and  the  results  of  the  dynamic  and 
wind  tunnel  tests  will  be  evaluated.  Included  in  this  evaluation  will  be  test 
to  analysis  comparisons  and  investigation  of  potential  regions  of  chaotic 
response  characteristics. 


5470D/6 


9 


As  discussed  in  Section  1.0,  the  nonlinear  system  response  of  an 
aerosurface  with  discrete  nonlinear  structural  elements  is  not  thoroughly 
understood.  The  objectives  of  the  program  plan.  Section  2.1,  is  to  perform 
research  to  develop  further  understanding  of  this  nonlinear  response 
phenomenon.  A  thirty-eight  month  analytical  and  experimental  research 
program,  including  documentation,  is  being  performed  to  meet  these 
objectives.  As  presented  in  Figure  6,  the  program  schedule  has  been 
established  to  provide  a  logical  flow  of  activities  through  the  program. 

MDMSC  has  responsibility  for  the  overall  program  with  UMR  serving  as  a 
subcontractor  to  MDMSC.  Roles  and  program  activities  for  UMR  and  MDMSC  are 
shown  in  Figure  7.  In  addition  to  overall  program  direction,  MDMSC  is 
performing  the  analytical  studies  and  assisting  with  test  data  analysis.  The 
design,  fabrication,  and  testing  of  the  aerosurface  models  is  being  performed 
by  UMR. 
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Substantial  progress  has  been  made  during  the  Basic  Year  towards 
achieving  the  overall  research  objectives.  Several  analysis  techniques  have 
been  established  and  used  to  obtain  further  understanding  of  the  nonlinear 
dynamic  response  characteristics  of  the  rigid  aerosurface  system.  Wind  tunnel 
testing  has  begun  for  the  rigid  aerosurface  and  a  majority  of  the  equipment 
necessary  for  dynamic  testing  has  been  acquired.  Comparison  of  analytical  and 
wind  tunnel  test  results  has  begun  for  the  rigid  aerosurface. 

3.1  Analysis 

A  variety  of  analysis  procedures  have  been  developed  during  the  Basic 
Year,  or  first  year  of  the  research  program.  These  are:  1)  a  flutter 
stability  analysis  utilizing  steady  aerodynamic  theory,  2)  a  transient 
response  analysis  using  steady  aerodynamic  theory,  and  3)  a  transient  response 
analysis  using  unsteady  aerodynamic  theory.  A  program  has  been  written  which 
computes  the  correlation  dimension  for  a  given  time  history.  In  addition,  the 
aerosurface  system  equilibrium  points  have  been  analytically  determined  for 
various  combinations  of  linear  and  nonlinear  root  support  springs. 

All  of  these  analysis  activities  are  directed  toward  the  objective  of 
understanding  the  nature  of  the  dynamic  response  of  a  aerosurface  containing 
structural  nonlinearities.  These  analysis  accomplishments  are  discussed  in 
detail  in  the  following  sections. 

Uol _ stability  Analysis  With  Steady  Aerodynamics  -  A  steady  aerodynamic 

flutter  analysis  was  performed  to  size  the  rigid  aerosurface  and  root  support 
mechanism  for  eventual  dynamic  and  wind  tunnel  testing.  The  sizing  was 
performed  to  assure  that  the  test  article  could  flutter  within  the  speed  range 
of  the  UMR  subsonic  wind  tunnel.  Analysis  results  defined  parametric 
relationships  between  rigid  aerosurface  weight,  inertia,  linear  support  spring 
stiffnesses,  and  location  of  the  elastic  axis,  aerodynamic  center,  and  center 
of  gravity.  Several  combinations  were  found  where  flutter  was  predicted  at 
speeds  in  the  100  to  150  ft/sec  range,  ideal  for  UMR  wind  tunnel  testing. 

These  results  were  used  as  guidelines  in  designing  the  test  article. 

The  steady  aerodynamic  flutter  method  was  chosen  for  use  in  the 
parametric  sizing  studies  because  the  method  has  proven  to  be  computationally 
efficient  and  reasonably  accurate  for  aerosurface  systems  similar  to  the  one 
of  this  study.  The  steady  aerodynamic  flutter  method  is  discussed  in 
References  6  through  8.  The  basic  assumption  made  in  this  method  is  that  the 
aerodynamic  lift  force  is  proportional  to  and  in  phase  with  the  angle  of 
attack.  Flutter  occurs  at  frequency  coalescence  for  the  two  degree  of  freedom 
representation  of  the  rigid  aerosurface.  This  leads  to  a  simple  closed  form 
definition  of  the  flutter  solution. 

The  system  equations  of  motion  for  the  rigid  aerosurface  in  Figure  5  are 
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The  two  degrees  of  freedom  are  aerosurface  root  pitch,  a,  and  plunge,  h, 
shown  in  Figure  5.  The  static  unbalance  Ta  is: 


Ta  =  mad  ( 2  ) 

For  low  speed,  incompressible  steady  aerodynamics,  the  aerodynamic  lift 
force  term  is 

L  =  qSCj^  a  (3) 

where  we  have  the  dynamic  pressure  and  lift-slope  coefficient  given  as 
q  =  1/2  pV2 

Cut  =  2ir  <4> 

The  aerodynamic  moment  is  given  as 

M  =  L( a  +  l/2b)  (5) 


Here  it  has  been  assumed  that  the  aerodynamic  center  is  at  the  quarter  chord. 
Figure  5. 

Substituting  the  aerodynamic  terms  into  the  system  equations.  Equation  (1) 
becomes 

mhh  +  Ta  a  +  Khh  +  1/2  pV^C^a  =  0  (6) 

Tah  +  Ia  a  ♦  [Ka  -  1/2  pV2SCLa(a  +  l/2b)]a  =  0 

Assuming  harmonic  motion  and  for  non-trivial  solutions,  the  determinant 
of  coefficients  of  Equation  (6)  must  be  zero.  Following  the  approach  of 
Reference  8,  and  defining  the  following  terms 

A  =  ir  pV2S 

3  =  nh1*  ~  Ta 
C  =  Ta  ♦  mh(a  +  l/2b) 

D  =  -Khla  -  Kamh  ( 7 ) 

E  =  -Kh  (a  +  l/2b) 

F  =  KhKa 

we  obtain  the  following  expression  for  the  parameter  A  at  flutter 


Two  values  for  Ap  are  obtained  with  the  lower  value  being  related  to  the 
system  flutter  velocity  by 


Example  results  from  the  steady  aerodynamic  flutter  solution  are  shown 
in  Figures  8  through  11.  For  these  cases,  the  aerosurface  was  modeled  as  a 
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rigid  plate.  The  center  of  gravity  was  held  constant  at  the  half-chord.  The 
total  weight  of  the  aerosurface  and  the  support  mechanism  was  assumed  to  be  5 
pounds.  The  aerodynamic  center  was  coincident  with  the  elastic  axis.  The 
aerodynamic  center  and  elastic  axis  move  from  the  quarter-chord  aft  towards 
the  half-chord.  The  parameter  xa  is  a  nondimensional  form  of  d.  Figure  5, 
such  that 

d  =  b  xa  (10) 

The  term  R  is  the  zero  airspeed  frequency  ratio  defined  by 

R  =  uh/ua  (11) 

Based  on  operational  characteristics  of  the  UMR  wind  tunnel,  the  goal 
for  wind  tunnel  testing  was  in  the  0.1  to  0.2  Mach  number  range.  From  the 
results  shown  in  Figures  8  through  11,  a  rigid  aerosurface  test  article 
weighing  about  5  pounds  with  pitch  and  plunge  frequencies  in  the  8  to  12  Hz 
range  seemed  ideal  for  purposes  of  this  study.  These  results  formed  the  basis 
of  the  eventual  definition  of  the  rigid  aerosurface  configuration. 

3.1.2  Transient  Analysis  With  Steady  Aerodynamics  -  A  steady  aerodynamic 
transient  response  simulation  has  been  developed.  The  program  has  been 
exercised  using  the  physical  properties  of  a  preliminary  aerosurface  test 
article  design.  The  program  will  be  used  in  the  upcoming  Option  I  year  to 
establish  conditions  for  dynamic  testing  and  to  obtain  results  for  comparison 
with  the  dynamic  test  results.  The  simulation  performs  transient  analysis  of 
the  two  degree  of  freedom  rigid  aerosurface.  Figure  5,  with  nonlinear  root 
supports  and  assuming  steady  aerodynamics.  The  simulation  uses  the  ACSL 
(Advanced  Continuous  Simulation  Language)  routine.  Reference  9. 

The  governing  second  order  equations  of  motion.  Equation  (6),  are 
transformed  to  first  order  equations  using  the  state  variable  transformation 
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The  system  of  first  order  equations  are  then  given  as 
M  X  +  K  X  =  0 
where  we  have 
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Note  that  velocity  proportional  damping  terms  and  Ca  has  been 
included  in  these  equations. 

The  ACSL  routine  is  used  to  numerically  integrate  Equation  (13).  A  flow 
velocity  is  assumed  and  the  system  is  given  a  prescribed  set  of  initial 
conditions.  The  resulting  time  solution  is  then  examined  for  stability  and 
chaos.  In  this  analysis,  the  root  springs  Kh  and  Ka  are  no  longer 
constants,  but  due  to  assumed  nonlinearities  are  dependent  on  aerosurface 
response.  During  a  simulation,  the  instantaneous  values  of  these  stiffnesses 
follow  the  relationships  shown  in  Figure  2. 

Due  their  relationship  with  the  dynamic  testing,  results  from  the  steady 
aerodynamic  transient  response  analysis  will  be  included  in  the  Option  I  year 
technical  report. 

3.1.3  Transient  Analysis  With  Unsteady  Aerodynamics  -  An  unsteady  aerodynamic 
transient  response  simulation  has  also  been  developed.  The  simulation  has 
been  used  to  establish  conditions  for  wind  tunnel  testing  and  to  obtain 
results  for  comparison  with  the  wind  tunnel  test  results.  The  simulation 
performs  transient  analysis  of  the  two  degree  of  freedom  rigid  aerosurface. 
Figure  5,  with  nonlinear  root  supports  and  assuming  unsteady  aerodynamics. 

The  simulation  uses  the  ACSL  (Advanced  Continuous  Simulation  Language)  routine. 

Two-dimensional,  incompressible  Theodorsen  aerodynamics.  References  10 
and  11,  are  used  for  this  transient  response  analysis.  Since  the  unsteady 
aerodynamic  coefficients  are  known  only  at  discrete  values  of  reduced 
frequency,  an  approximate  transfer  function  is  used  to  develop  the  continuous 
function  necessary  for  transient  analysis.  Additional  modifications  are 
incorporated  so  that  the  transient  flutter  solution  is  applicable  to  arbitrary 
motion  as  opposed  to  the  commonly  assumed  simple  harmonic  motion.  This  method 
is  described  in  detail  in  References  12  and  13. 

For  unsteady  aerodynamics,  the  equations  of  motion  for  the  rigid 
aerosurface.  Equation  (1),  are  expressed  as 

mh  h  ♦  Ta  a  +  Ch  h  ♦  Kh  h  =  L 

Ta  'h  +■  Ia  a  +  Ca  a  +  Ka  a  =  M  (15) 


where  the  Cxx  terms  are  the  unsteady  Th^^Hrir sen  coefficients.  The 
aerodynamics  are  approximated  by  Roger's  transfer  function.  Reference  14, 
which  consists  of  a  second  order  polynomial  in  terms  of  reduced  frequency  k 
with  a  finite  series  of  poles. 


[ C ( k ) ]  =  [A,]  +  [A  )  ( ik )  +  (A  1  ( ik ) 2  ♦ 
x  Z  3 


M  +  3 


[Al(ik) 

_ m _ 

(  ik )  +  k 


(  16) 


where  in  the  above  A^  is  the  noncirculatory  static  aerodynamics ,  A2(ik)  is 
( he  aerodynamic  damping,  Aj(ik)2  is  the  apparent  mass,  and  the  summation 
terms  are  circulatory  aerodynamic  lags. 
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The  value  of  m  can  be  1,  2,  3,  or  4,  depending  on  the  required 
accuracy.  All  analysis  performed  to  date  for  this  study  has  used  a  value  of  m 
equal  to  1.  The  transfer  function  circulatory  aerodynamic  lag  terms,  km’s, 
are  computed  from  a  best  fit  of  the  Theodorsen  function.  The  Aj^ 
coefficients  are  computed  by  the  least-squares  error  method  for  the  reduced 
frequencies  at  which  the  transfer  function  is  defined. 

For  the  aerosurface  system  in  this  study  we  have: 
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The  following  definitions  are  used. 


X  =  V 
•  •  • 
X  =  Y 


(CJ  =  [C]  -  q  V  [A2] 

[K]  =  [K]  -  q  [Aj]  (18) 

=  qP4  [A4] 

[Ax)  =  [A1 J  +  [A4J 


A  set  of  six  first  order  differential  equations  are  then  obtained  by 
utilizing  the  inverse  Laplace  transform,  the  response  to  a  unit  impulse,  and 
the  convolution  integral.  The  Z  terms  are  the  generalized  aerodynamic  lags. 
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The  system  equations.  Equation  (19)  have  been  programmed  in  ACSL  in  a 
manner  similar  to  that  for  steady  aerodynamics.  Aerosurface  transient 
response  is  obtained  for  varying  combinations  of  stiffness  nonlinearities, 
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initial  conditions,  and  flow  conditions.  As  with  the  steady  aerodynamic 
simulation,  Ka  and  are  variable  and  functions  of  the  aerosurface 
response.  Results  for  these  simulations  are  compared  to  wind  tunnel  test 
results  in  Section  3.3. 


3.1.4  Correlation  Dimension  -  An  automated  method  for  computing  the 
correlation  dimension  of  a  phase  space  trajectory  has  been  developed  to  aid  in 
identifying  chaotic  behavior.  The  correlation  dimension  is  a  probabilistic 
measure  which  can  be  used  to  classify  attractors  as  being  either  simple  or 
chaotic.  Simple  attractors  have  integer  dimensions  while  chaotic  attractors 
have  fract;onal  dimensions. 


The  correlation  dimension  formulation  used  in  this  study  is  taken  from 
Reference  15.  For  a  discretized  phase  space  trajectory  of  N  points,  a 
correlation  function  is  defined  as 


C(r)  = 


N 

r 


where  Nr  is  the  number  of  pairs  of  points  X^  and  Xj  such  that  |  X^  - 

Xj|  <  r.  This  function  exhibits  a  power  law  dependence  on  r  as  r  approaches 

zero  for  many  attractors  such  that 


Lim  C(r)  =  ar^  (21) 

r  — *-0 


for  some  constant  a.  The  correlation  dimension  dc  is  defined  using  the 
slope  of  the  log  C(r)  versus  log  r  curve 


d 

c 


log  r 


(22) 


The  correlation  dimension  of  the  Henon  map  has  been  computed  as  a  test 
case  for  these  calculations.  The  Henon  map  is  defined  as 

xn  +  l  =  x-0  -  axn  +  Yn  <23> 

Yn+1  =  b  xn 

For  this  test  case,  a  has  been  defined  as  1.4  and  b  as  0.3. 

A  total  of  5000  points,  N,  in  the  x-y  phase  space  were  computed  and  are 
shown  in  Figure  12.  The  log  C(r)  versus  log  r  plot  is  shown  in  Figure  13. 
Clearly  there  is  a  limited  range  where  the  slope  is  approximately  constant  and 
only  this  range  is  used  when  computing  the  correlation  dimension.  Computing 
the  slope  via  a  least  squares  fit  for  values  of  log  r  between  52  and  64  yields 
a  correlation  dimension  of  1.22.  This  compares  to  the  value  established  in 
Reference  16  of  1.21  ±  0.01. 

The  usefulness  of  the  correlation  dimension  in  classifying  attractors 
for  the  aerosurface  system  being  used  in  this  study  will  be  investigated 
during  Option  I  year.  The  current  plan  is  to  use  the  embedding  technique, 
also  known  as  the  pseudo-phase  space  approach,  as  described  in  Reference  17. 
This  technique  has  proven  to  work  well  with  experimental  data.  Reference  18. 
Instead  of  measuring  Xj^,  only  a  single  component  X^  is  measured.  A  new 
n-dimensional  phase  space  is  constructed  as: 

22 
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-i  ~  (Xj^,  X|  +  (-  X^+2t'  •••'  ^i-t-nt^ 

where  t  is  an  almost  arbitrarily  chosen  time  constant.  The  main  restriction 
on  choosing  t  is  to  avoid  the  system's  basic  oscillation  period. 


The  correlation  dimension  is  computed  using  the  Z^  vectors.  The  value 
of  the  embedding  dimension,  n,  is  increased  until  a  constant  value  of  the 
correlation  dimension  is  obtained.  This  procedure  will  be  used  to  assess  the 
possibility  of  chaos  in  both  the  analytical  transient  response  results  and  the 
experimental  data. 

3.1.5  Equilibrium  Point  Identification  -  The  equilibrium  points  for  systems 
with  various  nonlinear  root  spring  combinations  and  aerodynamic  center 
locations  have  been  analytically  determined.  Chaotic  behavior  of  the 
aerosurface  system  would  seem  to  be  caused  by  the  response  jumping  from  one 
equilibrium  point  to  another.  Several  nonlinear  root  spring  and  aerodynamic 
center  location  combinations  yielding  multiple  system  equilibrium  points  were 
found.  The  analytical  and  experimental  search  for  chaotic  response  will  be 
focused  on  these  combinations. 

Starting  with  the  equations  of  motion  for  the  aerosurface,  such  as 
Equation  (1),  the  dynamic  terms  are  set  equal  to  zero  to  determine  the 
equilibrium  points.  Thus  we  are  looking  for  solutions  to  the  following 
relationships 


(25) 

where  Fsh  and  Msa  are  the  forces  in  the  root  plunge  and  pitch  springs. 

To  simplify  the  analysis,  steady  aerodynamics  are  assumed.  The  lift  force 
at  the  elastic  axis  is  given  as 

L  =  1  s  cLa  a  (26) 

where 

q  =  1/2  p  (a0M)2  (27) 

cLa  =  2ir 


Fsh  (h)  =  -L(a) 
Msa  (a)  =  M(a) 


I 


Combining  Equations  (26)  and  (27),  we  define  A  as  the  aerodynamic  constant 

A  =  itSpa2  (28) 

The  aerodynamic  moment  at  the  elastic  axis  is  thus  expressed  as: 

M  =  Lr  =  ArM2a  (29) 


where  r  is  the  distance  from  the 
(positive  aft).  To  simplify  the 
considered.  This  corresponds  to 
with  the  aerodynamic  center  (r  = 
aerodynamic  center  (r  >  0).  The 


aerodynamic  center  to  the  elastic  axis 
analysis,  no  negative  r  values  are 
a  system  where  the  elastic  axis  is  coincident 
0)  or  the  elastic  axis  is  behind  the 
system  equations  become 


Fsh(h)  =  -AM2<x 

Msa  <a)  =  ArM2a  (30) 
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The  nonlinear  spring  force  and  moment  terms,  Fsj1(h)  and  Msa(a), 
are  now  considered.  For  the  nonlinear  springs  used  in  this  study,  the  load 
displacement  relationship  can  be  broken  into  linear  segments  as  shown  in 
Figure  2.  For  the  linear  spring  we  have 

L  =  kx  for  -“  £  x  £  00 

For  the  freeplay  spring,  we  have 

L  =  k(x  +  S)  for 

L  =  0  for 

L  =  k ( x  -  S)  for 

and  for  the  preload  spring,  this  becomes 

L  =  kx  for 

L  =  kP  for 

L  =  k ( x  -  2S )  for 

In  order  to  determine  the  equilibrium  points,  all  possible  displacement 
combinations  are  considered  for  a  given  pair  of  springs.  For  example,  a 
system  with  two  freeplay  springs  has  the  following  possible  displacement 
combinations . 
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There  are  three  possibilities  for  each  displacement  combination:  1)  a  unique 
equilibrium  point  exists  (the  system  is  determinate),  2)  an  infinite  number  of 
equilibrium  points  exist  (the  system  is  indeterminate),  or  3)  no  equilibrium 
point  exists  (the  system  is  undefined). 

A  second  assumption  used  in  this  analysis  was  that  only  air  stream 
velocities  below  the  system  divergence  velocity  are  considered.  Divergence  is 
a  special  case  of  an  equilibrium  point  where  the  displacement  is  equal  to 
infinity.  At  divergence,  the  aerodynamic  moment  equals  the  spring  moment. 

That  is 


£  x  <_  -S 
-S  £  x  <  S 
S  <  x  ~  - 


-00  £  X  £  P 
P  £  X  £  P  +  2S 
P  +  2S  <  x  <  00 


(31) 


(32) 


(33) 


Ka  =  ArM2  (34) 

The  equilibrium  points  for  several  spring  combinations  are  summarized  in 
Figure  14.  The  details  of  establishing  these  characteristics  are  presented  in 
Appendix  A. 


5547D/8 


26 


SUPPORT  TYPE 
2  LINEAR/NO  AERO  MOMENT 
2  LINEAR/WITH  AERO  MOMENT 
2  PRELOAD /NO  AERO  MOMENT 
2  PRELOAD/WITH  AERO  MOMENT 
2  FREEPLAY/NO  AERO  MOMENT 
2  FREEPLAY/WITH  AERO  MOMENT 


EQUILIBRIUM  POINTS 

DETERMINATE 

INDETERMINATE 

UNDEFINED 

1 

0 

0 

1 

0 

0 

1 

0 

8 

3 

0 

6 

2 

3 

4 

4 

1 

4 
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3.2  Experimentation 

The  majority  of  test  activity  during  the  Basic  Year  dealt  with  wind  tunnel 
testing  of  the  rigid  aerosurface.  The  emphasis  of  this  report  will  be  on  this 
testing.  For  completeness,  a  summary  is  given  of  the  design  and  initial 
fabrication  of  the  rigid  aerosurface  dynamic  testing  equipment.  The  dynamic 
testing  will  be  performed  during  the  upcoming  Option  I  year. 

3.2.1  Wind  Tunnel  Testing  -  The  overall  wind  tunnel  test  facility  is  shown  in 
Figure  15.  The  UMR  subsonic  wind  tunnel  is  a  closed  circuit,  single  return 
tunnel.  The  test  section  has  an  internal  width  of  4  feet,  height  of  3.67 
feet,  and  a  length  of  11  feet.  It  is  constructed  of  plywood  with  Plexiglas 
sides.  The  tunnel  has  a  settling  chamber  with  a  counteraction  ratio  of  9  to  1 
and  two  screens  to  equalize  the  tunnel  velocity  profile  and  reduce  turbulence 
intensity.  Both  screens  are  a  square  14  mesh  (14  wires  per  inch  across  the 
mesh)  with  a  wire  diameter  of  0.020  inch  for  the  upstream  screen  and  0.022 
inch  for  the  downstream  screen.  Open  area  of  the  screen  is  about  70%  of  total 
area.  Both  the  open  area  ratio  of  the  screens  and  the  counteraction  ratio 
follow  Bradshaw's  recommendations  for  boundary  layer  testing. 

The  power  plant  for  the  wind  tunnel  is  a  350  horsepower  diesel  engine. 

This  engine  drives  a  variable  displacement  hydraulic  pump  which  in  turn  drives 
a  fixed  displacement  hydraulic  motor.  The  motor  is  directly  connected  to  the 
wind  tunnel  fan  through  a  driveshaft  provided  with  universal- joint  couplings 
at  each  end.  The  airflow  power  system,  through  its  variable  pitch  fan  blade 
settings  and  variable  speed  hydraulic  transmission,  provides  a  wide  range  of 
tunnel  velocities  with  variable  control  of  velocity  within  a  range  allowed  by 
the  fan  pitch.  Feedback  command  control  of  fan  speed  is  performed  by  a  Moog 
controller  and  a  D.C.  tachometer.  A  speed  range  from  near  zero  to  in  excess 
of  200  miles  per  hour  is  possible. 

As  seen  in  Figure  15,  a  support  frame  is  attached  to  the  floor  and 
extended  around  the  outside  of  the  tunnel.  The  aerosurface  is  supported  at 
the  top  of  the  wind  tunnel  by  a  root  support  mechanism  mounted  to  the  support 
frame.  Additional  details  of  the  aerosurface  installation  are  shown  in 
Figures  16  and  17.  The  aerosurface  support  shaft.  Figure  16,  extends  through 
the  tunnel  ceiling  to  the  support  mechanism  shown  in  Figure  17.  The  suppoit 
mechanism  is  free  to  move  in  the  horizontal  direction  on  four  linear  bearings 
supported  by  two  precision  shafts.  This  horizontal  motion,  aerosurface  plunge 
degree  of  freedom  h,  is  monitored  by  a  LVDT  as  shown  in  Figure  18.  Rotational 
motion  of  the  aerosurface  support  shaft  corresponds  to  the  pitch  degree  of 
freedom  a.  A  RVDT,  Figure  18,  monitors  this  degree  of  freedom. 

I  Output  from  these  transduces  go  to  signal  conditioners  and  then  to  a 

digital  processing  oscilloscope.  Figure  19.  With  this  oscilloscope,  the 
digitized  response  data  may  be  displayed  as  time  history  information.  In 
addition,  phase  planes  and/or  frequency  domain  information  may  be  obtained  for 
these  time  history  data.  The  digitized  data  may  be  down-loaded  to  the  data 
storage  personal  computer  for  paper  plotting  and/or  further  analysis. 

The  aerosurfaces  root  stiffness  and  associated  nonlinearities  may  be 
varied  thro"gh  features  incorporated  in  the  design  of  the  support  mechanism. 
Figure  20.  Steel  leaf  spring  material  of  varying  thicknesses  are  used  for  the 
pitch,  Kh,  and  plunge,  Ka,  springs.  Various  magnitudes  of  freeplay  gaps 
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°°ilUPPORT  mechanism  . 
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FIGURE  16.  AEROSURFACE 
INSTALLED  IN  WIND  TUNNEL 
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FIGURE  17.  ROOT  SUPPORT  MECHANISM 


FIGURE  18.  INSTRUMENTATION 


FIGURE  19.  TEST  CONTROL  AND  DATA  ACQUISITION  SYSTEM 
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are  obtained  by  adjusting  the  pitch  and/or  plunge  support  stops  shown  in 
Figure  20.  A  linear  stiffness  is  obtained  by  tightly  securing  the  stops 
|  against  the  leaf  spring. 

I 

Mass  and  inertia  properties  of  the  moving  portions  of  the  aerosurface  and 
1  support  mechanism  have  been  determined.  Several  different  combinations  of 

spring  elements  are  available  to  obtain  various  combinations  of  root  support 
stiffnesses.  A  summary  of  this  mass  and  stiffness  information,  in  addition  to 
^  other  aerosurface  data,  is  given  in  Figure  21. 

i  3.2.2  Dynamic  Testing  -  Dynamic  testing  is  planned  during  the  Option  I  year 

to  obtain  additional  data  to  complement  the  wind  tunnel  data  for  the  nonlinear 
,  system.  The  dynamic  testing  will  employ  an  equivalent,  mechanically  excited 

:  system.  During  this  testing,  excitation  will  be  controlled  by  a  "feed  back" 

system,  in  a  manner  analogous  to  the  aerodynamic  forcing  function  associated 
|  with  aerosurface  motion  in  a  subsonic  airstream. 

I  In  the  experimental  approach  for  dynamic  testing.  Figure  22,  a  transducer 

will  measure  the  root  pitch  parameter  a  and  feed  this  information  back  as 
the  controlling  function  for  a  mechanical  actuator.  The  signal  from  the  pitch 
degree  of  freedom  RVDT  goes  to  a  signal  conditioner  and  then  to  the  actuator 
controller.  Also  fed  to  the  controller  is  the  output  from  a  load  cell  on  the 
actuator  drive  shaft.  The  controller  maintains  an  actuator  load  which  is 
equal  to  the  aerodynamic  lift  associated  with  the  pitch  response  for  an 
assumed  flow  condition. 

A  hydraulic  actuator  configuration  has  been  selected  for  this  testing. 
Thus,  the  controller  governs  a  servovalue  which  controls  the  hydraulic  flow  to 
a  piston.  A  majority  of  individual  items  have  been  assembled  for  this  test 
setup.  The  hydraulic  servovalve  is  on  order  and  the  actuator  piston  is  being 
procured  from  an  outside  machine  shop.  There  has  been  a  delay  in  the  delivery 
of  these  items.  Therefore  the  rigid  aerosurface  dynamic  testing  will  be 
performed  later  in  Option  I  year  following  the  wind  tunnel  testing. 

_ Test  to  Analysis  Comparison 

Wind  tunnel  testing  of  the  rigid  aerosurface  system  has  begun  with  several 
preliminary  cases  having  been  run.  A  number  of  these  cases  are  described  in 
this  section.  Time  histories  of  the  experimental  data  are  shown.  The 
experimental  results  are  compared  with  analytical  results  from  the  transient 
response  simulations  using  unsteady  aerodynamics  as  described  in  Section 
3.1.3.  The  preliminary  test  results  will  be  used  to  determine  the  matrix  of 
conditions  for  the  rigid  aerosurface  complete  test  program  to  be  completed 
during  the  Option  I  year.  Lessons  learned  while  performing  these  initial 
tests  will  be  used  to  refine  the  test  procedure  prior  to  performing  the 
complete  test  program. 

All  of  the  test  results  described  in  this  section  have  a  plunge  spring 
constant  of  37.8  lb/in  and  a  pitch  spring  constant  of  704.4  in-lb/rad.  These 
values  correspond  to  a  plunge  frequency  of  8.21  Hz  and  a  pitch  frequency  of 
10.2  Hz. 
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ProDertv 

Symbol 

Value 

Semichord  length 

b 

3.5  in 

Span 

- 

24.5  in 

Distance  from  midchord  to 
elastic  axis  (+  aft) 

a 

-1.75  in 

Distance  from  elastic  axis  to 
center  of  gravity  (+  aft) 

d 

0.84  in 

Pitch  mass 

ma 

5.037(10~3)  lb-sec2/in 

Plunge  mass 

mh 

14 . 22 ( 10-3 )  lb-sec2/in 

Pitch  moment  of  inertia 
about  elastic  axis 

la 

0.1715  lb-sec2  -  in 

Pitch  linear  spring  constant 
range 

ka 

732-1831  in-lb/rad 

Plunge  linear  spring  constant 
range 

kh 

13.7-215  lb/in 

Pitch  spring  deadband  (max) 

+Sa 

11  degrees 

Plunge  spring  deadband  (max) 

±Sh 

0.56  in 

FIGURE  21.  RIGID  AEROSURFACE  PROPERTIES 
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(a)  IDEALIZATION 


(b)  IMPLEMENTATION 

FIGURE  22.  DYNAMIC  TEST  IMPLEMENTATION 
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The  first  case  is  for  a  linear  system.  Time  histories  of  *-he  experimental 
data  are  shown  in  Figure  23.  This  data  is  for  a  velocity  of  122  ft/sec.  The 
system  is  close  to  flutter  at  this  velocity,  as  can  be  seen  from  the 
relatively  low  damping  in  the  experimental  time  histories.  Additional 
experimentation  has  determined  the  flutter  speed  to  be  between  125  ft/sec  and 
130  ft/sec.  Analytical  results  for  this  case  are  shown  in  Figures  24  and  25. 
From  these  results  it  is  evident  that  analytical  flutter  occurs  between  110 
and  115  ft/sec.  Work  is  in  progress  to  determine  why  the  analytically 
determined  flutter  speed  is  lower  than  the  flutter  speed  obtained  from 
testing.  One  possible  reason  is  that  zero  mechanical  damping  was  included  in 
the  analysis.  An  appropriate  value  for  mechanical  damping  will  be  determined 
and  incorporated  in  the  analysis  during  subsequent  analyses. 

The  second  set  of  experimental  data  is  for  a  system  which  has  a  linear 
plunge  spring  and  a  +2.5  degree  freeplay  deadband  in  the  pitch  degree  of 
freedom.  Time  histories  for  the  experimental  data  are  shown  in  Figure  26. 

The  system  is  experiencing  a  stable  limit  cycle  oscillation  of  approximately 
+15  degrees  pitch  and  +1  inch  plunge  motion.  These  data  are  shown  at  a  time 
following  the  startup  transient  in  the  test  system.  The  analytical  time 
histories  are  shown  in  Figure  27.  After  a  startup  transient  the  analysis  also 
predicts  stable  limit  cycle  response.  In  this  case  the  analysis  predicts 
limit  cycle  oscillation  of  approximately  +23  degrees  pitch  and  +1.6  inches 
plunge.  Reasons  for  the  difference  between  the  amplitudes  of  limit  cycle 
oscillations  for  the  experimental  and  analytical  data  are  being  investigated. 

A  possibility  is  that  as  in  the  linear  case,  zero  mechanical  damping  was  used 
in  the  analysis.  This  will  be  studied  during  the  Option  I  year.  To  aid  in 
understanding  the  relation  between  experimental  and  analysis  results,  the 
frequency  content  of  the  data  will  also  be  compared. 

The  third  set  of  data  is  for  a  system  with  a  +2.5  degree  freeplay  deadband 
in  pitch  and  a  ±0.25  inch  freeplay  deadband  in  plunge.  Experimentally 
obtained  time  histories  are  shown  in  Figures  28  and  29.  At  the  lower  velocity 
of  61  ft/sec,  the  system  is  damped  and  the  response  dies  out.  A  stable  limit 
cycle  oscillation  is  obtained  for  the  higher  velocity  of  81  ft/sec.  The 
analytical  time  histories  for  the  81  ft/sec  case  are  shown  in  Figure  30.  As 
can  be  seen  the  nature  of  the  experimental  and  analytical  response  is 
different.  In  particular,  the  analytical  response  seems  to  be  beating  in  both 
degrees  of  freedom  with  a  slight  indication  of  beating  in  the  experimental 
pitch  degree  of  freedom.  The  cause  of  this  phenomena  is  not  understood  and 
will  be  investigated  during  the  next  year. 
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4.0  OPTION  I  YEAR  PLANS 


During  the  upcoming  Option  I  year,  this  research  program  will  build  on  the 
accomplishments  of  the  completed  Basic  Year.  A  schedule  of  these  activities 
is  presented  in  Figure  6,  Section  2.2. 

Extensive  wind  tunnel  testing  will  be  conducted  for  the  rigid  aerosurface 
with  varying  combinations  of  root  support  nonlinearities.  Aerosurface 
response  characteristics  will  be  determined  as  a  function  of  air  speed, 
magnitude  and  type  of  nonlinearities,  and  initial  conditions.  Post-test 
analysis  of  measured  response  will  focus  on  identifying  potential  regions  of 
chaotic  response.  Included  will  be  frequency  and  power  spectral  density 
analyses,  plotting  of  phase  planes  and  Poincare’  maps,  and  determination  of 
the  associated  correlation  dimension. 

Numerical  simulations  of  the  nonlinear  dynamics  problem  will  be  conducted 
in  parallel  with  the  wind  tunnel  testing.  These  simulations  will  employ  the 
unsteady  aerodynamic  representation.  Section  3.1.3.  Aerosurface  response 
characteristics  will  be  evaluated  numerically  to  determine  the  nature  of  this 
response  in  terms  of  initial  conditions,  air  speed,  and  nonlinearities. 

Results  from  these  numerical  studies  will  aid  in  defining  wind  tunnel  test 
conditions.  In  addition,  these  simulation  results  will  be  evaluated  for  the 
potential  of  chaotic  behavior  in  a  maneuver  similar  to  that  used  for  the  wind 
tunnel  test  data. 

In  reviewing  results  from  the  wind  tunnel  testing,  it  became  apparent  that 
there  was  slippage  in  the  root  spring  of  the  pitch  degree  of  freedom.  For  the 
current  design,  this  leaf  spring.  Figure  20,  is  held  in  place  with  set 
screws.  During  testing  there  appeared  to  be  relative  motion  through  this 
attachment.  Early  in  the  Option  I  year  the  attachment  of  the  pitch  spring 
will  be  modified  to  ensure  a  positive  support  at  this  point. 

During  the  coming  year,  the  dynamic  test  setup  will  be  completed  and 
dynamic  testing  begun  for  the  rigid  aerosurface.  This  testing  is  to  obtain 
data  to  augment  the  information  obtained  during  wind  tunnel  testing.  This 
information  will  enhance  the  understanding  of  the  behavior  of  the  nonlinear 
aeroelastic  system.  In  conjunction  with  this  testing,  numerical  simulation 
studies  will  be  performed  using  the  techniques  discussed  in  Section  3.1.2. 
These  analyses  will  employ  a  steady  aerodynamics  representation  which  is 
consistent  with  the  planned  dynamic  test  techniques  discussed  in.  Section 
3.2.2. 


UMR  has  completed  an  initial  design  and  assembled  a  prototype  model  for 
the  flexible  aerosurface.  This  is  well  ahead  of  the  scheduled  milestones 
shown  in  Figure  6.  The  initial  design/fabrication  activities  were  the  result 
of  a  senior  design  project  under  the  direction  of  Dr.  Eversman.  At  this  time. 
Dr.  Galecki  has  begun  fabrication  of  an  improved  prototype  of  the  flexible 
aerosurface  spar.  In  all  likelyhood,  final  fabrication  and  initiation  of  wind 
tunnel  testing  for  the  flexible  aerosurface  will  be  accomplished  during  the 
Option  I  year.  The  probability  of  this  activity  occuring  during  Option  I  has 
been  reflected  in  the  Figure  6  Program  Schedule. 
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5.0  CONCLUSIONS 


Substantial  progress  has  been  made  during  the  Basic  Year  towards  achieving 
the  overall  research  objectives.  These  objectives  are  to  obtain  an  increased 
understanding  of  the  dynamic  response  phenomenon  for  aerosurfaces  with 
discrete  structural  nonlinearities  within  the  context  of  chaotic  behavior. 

Several  analysis  techniques  have  been  established  and  used  to  obtain 
further  understanding  of  the  nonlinear  dynamic  response  characteristics  of  the 
rigid  aerosurface  system.  Wind  tunnel  testing  has  begun  for  the  rigid 
aerosurface  and  a  majority  of  the  equipment  necessary  for  dynamic  testing  has 
been  acquired.  Comparison  of  analytical  and  wind  tunnel  test  results  has 
begun  for  the  rigid  aerosurface. 

For  the  rigid  aerosurface,  similar  type  motion;  damped  decay,  stable  limit 
cycle,  and  flutter;  has  been  observed  in  both  the  analysis  and  experimental 
results.  A  number  of  response  characteristics  seen  in  the  analysis  results 
need  experimental  confirmation.  These  include  beat  phenomenon,  dependance  on 
initial  conditions,  and  indications  of  chaotic  response.  In  addition,  the 
influence  of  mechanical  damping  on  analysis  results  requires  further 
investigation.  These  topics  will  be  addressed  during  the  coming  year. 
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A  summary  of  the  professional  personnel  associated  with  the  research 
program  is  presented  in  the  following  section.  This  is  followed  by  a 
discussion  of  interactions  with  other  organizations. 

6.1  Personnel 

Dr.  Robert  M.  Laurenson  of  MDMSC  is  Principal  Investigator  and  is 
responsible  for  overall  program  management  and  performance  including  cost  and 
schedule  control.  Mr.  Anthony  J.  Hauenstein  is  Co-Investigator  reporting 
directly  to  Dr.  Laurenson  and  is  responsible  for  MDMSC  program  technical 
performance.  Dr.  Waiter  Eversman  is  Co-Principal  Investigator  and  directs  all 
program  activities  accomplished  at  UMR .  Assisting  Dr.  Eversman  at  UMR  are 
Dr.  Grzegorz  Galecki,  Dr.  Richard  T.  Johnson,  and  Mr.  Iyad  Khalil  Qumei. 

During  the  research  project  Dr.  Laurenson  reports  to  Dr.  John  L.  Gubser 
who  is  head  of  the  Structures  Department  at  MDMSC.  Due  to  his  background  as 
Principal  Investigator  during  the  previous  two  AFOSR  studies.  References  1  and 
2,  Dr.  Gubser  serves  as  a  Consultant  during  this  study.  In  this  capacity  he 
participates  in  program  reviews. 

Summary  information  for  these  individuals  follows. 


ROBERT  M.  LAURENSON 
Principal  Investigator 


Position:  Branch  Chief 

Structures  Department 

McDonnell  Douglas  Missile  Systems  Company  (MDMSC) 

St.  Louis,  Missouri 

Education:  1961;  B.S.,  Mechanical  Engineering,  University  of  Missouri-Rolla, 

Rolla,  MO 

1962;  M.S.,  Mechanical  Engineering,  University  of  Michigan, 

Ann  Arbor,  MI 

1969;  Ph.D.,  Mechanical  Engineering,  Georgia  Institute  of 
Technology  Atlanta,  GA 

Thesis  Title:  "Dynamic  Characteristics  of  Nonlinearity 
Damped  System" 


ANTHONY  J.  HAUENSTEIN 
Co -Investigator 


Position:  Engineer 

Structures  Department 

McDonnell  Douglas  Missile  Systems  Company  (MDMSC) 

St.  Louis,  Missouri 

Education:  1984;  B.S.,  Aerospace  Engineering,  Iowa  State  University,  Ames,  IA 

I 


5470D/11 


WALTER  EVERSMAN 
Co-Principal  Investigator 


Position: 


Education : 


Position: 


Education : 


Position: 


Education : 


Curators'  Professor 

Mechanical  and  Aerospace  Engineering  and  Engineering 
Mechanics  Department 

University  of  Missouri-Rolla  (UMR) 

Rolla,  Missouri 

1959;  B.S.,  Aerospace  Engineering;  Purdue  University, 

West  Lafayette,  IN 

1961;  M.S.,  Engineering  Mechanics;  Stanford  University, 
Stanford,  CA 

1964;  Ph.D.,  Aeronautics  &  Astronautics,  Stanford  University 
Stanford,  CA 

JOHN  L.  GUBSER 
Program  Consultant 

Manager 

Structures  Department 

McDonnell  Douglas  Missile  Systems  Company  (MDMSC) 

St.  Louis,  Missouri 

1958;  B.S.,  Mechanical  Engineering;  Washington  University, 

St.  Louis,  MO 

1959;  M.S.,  Engineering  Mechanics;  University  of  Illinois, 
Urbana,  IL 

1961;  Ph.D.,  Engineering  Mechanics;  University  of  Illinois, 
Urbana,  IL 

1961;  Fullbright  Fellowship,  Technological  University,  Delft  in 
the  Netherlands 


RICHARD  T.  JOHNSON 
Test  Instrumentation 

Professor 

Mechanical  and  Aerospace  Engineering  and  Engineering  Mechanics 
Department 

University  of  Missouri-Rolla  (UMR) 

Rolla,  Missouri 

1962;  B.S.,  Mechanical  Engineering;  University  of  Missouri  Rolla, 
Rolla,  MO 

1964;  M.S.,  Mechanical  Engineering;  University  of  Missouri-Rolla, 
Rolla,  MO 

Thesis  Title:  "The  Design,  Construction  and  Operation 
of  a  Conducting  Field  Analog  Computer" 

1968;  Ph.D.,  Mechanical  Engineering;  University  of  Iowa, 

Iowa  City,  IA 

Thesis  Title:  "Fluid  Jet  Transmission  of  a  Pressure 
Signal  from  a  Rotating  Shaft  to  a  Fixed  Receiver" 
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GRZEGORZ  GALECKI 
Mechanical  Hardware  Fabrication 


Position:  Visiting  Scholar 

Rock  Mechanics  &  Explosives  Research  Center 
High  Pressure  Waterjet  Laboratory 
University  of  Missouri-Rolla 
Rolla,  Missouri 

Education:  1974;  B.S.,  Master  of  Mechanical  Engineering,  Wroclaw  Technical 

University,  Wroclaw,  Poland 

1978;  Ph.D.  in  Technical  Sciences,  Wroclaw  Technical  University, 
Wroclaw,  Poland 

Thesis  Title:  "Mechanism  of  Cutting  of  a  Material  by 
High  Speed  Water  Jet” 


IYAD  KHALIL  QUMEI 

Coordination  of  Hardware  Fabrication  and  Testing 
Position:  Graduate  Student 

Mechanical  and  Aerospace  Engineering  Mechanics  Department 
University  of  Missouri-Rolla 
Rolla,  Missouri 

Education:  1986;  B.S.,  Mechanical  Engineering,  Yarmouk  University, 

Irbid,  Jordan 

1986;  Graduate  Student,  Mechanical  and  Aerospace  Engineering, 
University  of  Missouri-Rolla,  Rolla,  MO 

6.2  Interactions 

During  the  performance  of  this  research  program,  industry/university 
information  exchange  has  developed  in  addition  to  that  required  by  the 
research  objectives.  This  experience  has  been  selected  for  incorporation  in 
the  proceedings  of  the  ASME  1989  Mechanical  Engineering  Department  Heads 
Conference  which  was  held  March  8-10,  1989  in  Orlando,  Florida.  A  description 
of  this  interaction  is  contained  in  Appendix  B. 

Dr.  Eversman  has  obtained  $25K  in  research  funding  from  the  state  of 
Missouri  under  the  Missouri  Research  Assistance  Act  (MRAA) .  These  monies  are 
available  for  faculty  members  at  Missouri  institutions  who  are  doing  research 
for  Missouri  firms.  This  funding  is  a  direct  augmentation  of  the  AFOSR  funds 
that  UMR  is  receiving  via  their  subcontract  with  MDMSC. 

Briefings  have  been  made  to  a  number  of  organizations  which  included  a 
review  of  the  work  being  done  in  this  research  program.  A  listing  of  the 
briefings  held  during  the  Basic  Year  is  presented  in  Figure  31.  A  review  of 
the  research  program  objectives/results  was  recently  held  with  members  of  the 
Structural  Dynamics  and  Loads  Department  from  McDonnell  Aircraft  Company.  In 
addition,  periodic  status  briefings  have  been  made  to  MDMSC  management 
including  the  Vice  President  of  Engineering  and  the  Director  of  Design  and 
Technology  Engineering. 
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DATE 


INDIVIDUAL* S) 


ORGANIZATION 


14  April  1988 


12  July  1988 


17  August  1988 


31  August  1988 
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Dr.  Gerald  Seidel 
J.V.  Marron 
Jay  Schneider 
Michelle  Finegan 
Alkis  Koutsiouroubas 
Robert  Hay 


Structural  Design  &  Analysis  Branch 
Code  6041 

Naval  Air  Development  Center 
Warminister,  PA  18974 


Dr.  Daniel  R.  Mulville  Office  of  Aeronautic  &  Space 

Technology 

National  Aeronautics  &  Space 
Administration 
Washington,  D.C.  20546 


Randall  C.  Davis  Thermal  Structures  Branch 

John  Shideler  NASA  Langley  Research  Center 

H.  Neale  Kelly  (PRC)  Hampton,  VA  23665 

Wayne  Sawyer 

Charles  Camarda 

Jim  Robinson 

Vic  Spain  (PRC) 

Jeff  Cerro  (PRC) 

Stuart  Jones  (PRC) 

Shawn  Birchenough 
Kevin  Rivers 
Steve  Scotti 
Lynn  Bownan 
David  Glass 


Dr.  Craig  Porter 
Michael  Munson 
Jim  McManigal 
Ted  Hicks 
Linda  Finco 


Thermal/Structures  Branch 
Code  3242 

Naval  Weapons  Center 
China  Lake,  CA  93555 


FIGURE  31.  BASIC  YEAR  SUMMARY  BRIEFINGS 
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EQUILIBRIUM  POINT  IDENTIFICATION 
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This  Appendix  contains  examples  which  demonstrate  how  the  equilibrium 
points  discussed  in  Section  3.1.5  were  identified.  All  the  displacement 
combinations  for  the  two  cases  of:  1)  two  linear  springs  with  no  aerodynamic 
moment  and  2)  two  freeplay  springs  with  an  aerodynamic  moment  are  shown. 
Equilibrium  points  for  systems  with  other  spring  combinations  may  be  defined 
in  a  similar  manner.  Figure  2  from  Section  1.0  is  repeated  here  as  Figure  A-l. 

Case  1  -  Two  Linear  Springs  With  No  Aerodynamic  Moment 

This  case  has  one  possible  set  of  displacement  combinations  defined  as 

-  00  <  a  <  00  and  -  00  <  h  <  115 

For  this  case  the  system  equations,  obtained  from  Equation  (30),  are 

Kh  =  -AM2a 
Ka  a  =  0 

Solving  the  above  for  h  and  a  yields: 

a  =  0  and  h  =  0 


This  is  a  unique  equilibrium  point  and  the  system  is  determinate. 

CLa.se  .2  -  Two  Freeplay  Sptisgs  with  An  .Aerodynamic  Moment 

This  case  has  nine  possible  displacement  combinations.  Each  of  these 
combinations  are  individually  considered. 


-  00  <  <x  <  -Sa  and  -  00  <  h  <  -S^ 

The  system  equations  for  this  situation  are 

Kh  h  +  Rh  Sh  =  -AM2a 
Ka  a  +  Ka  Sa  =  AM2ra 

Solving  for  h  and  a  yields 

a  =  -Ka  Sa/(Ka  -  AM2r) 
h  =  -AM2a/Kh  -  Sh 

From  the  displacement  assumption,  a  is  always  negative.  This  causes  h  to 
always  be  greater  than  S^,  which  violates  the  displacement  assumption  for 
h.  Thus,  Combination  1  has  no  equilibrium  point  and  the  system  is  undefined. 

o  Combination  2: 

-Sa  <  a  <  Sa  and  -  00  <  h  <  -S^ 

The  system  equations  are 
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Solving  for  h  and  a  yield 


a  =  0  and  h  =  -Sa 


This  is  a  unique  equilibrium  point  and  the  system  is  determinate, 
o  Combination  3; 

Sa  <  a  <  00  and  -  ®  <  h  <  -Sh 
The  system  equations  are 


Khh  +  KhSh  =  -AM2a 
Ka  a  -  KaSa  =  AM2  rd 

Solving  for  h  and  yields 

a  =  K^/O^  -  AM2r) 
h  =  -AM2a/Kh  -  Sh 


This  is  a  unique  equilibrium  point  and  the  system  is  determinate, 
o  Combination  4: 

-  00  <  a  <  -Sa  and  -Sh  <  h  <  Sh 
The  system  equations  are 
0  =  -AM2a 

Kaa  +  KaSa  =  AM2ra 

These  equations  yield  two  solutions  for  a  which  are 


a  =  0 

a  =  -KaSa/(Ka  -  AM2r) 

and  thus  this  combination  has  no  solution  and  the  system  is  undefined, 
o  Combination  5; 

-  sa  < a  <  sa  and  _sh  <  h  <  Sh 

The  system  equations  are 

0  =  -AM2a 
0  =  AM2  ra 


Solving  for  a  yields 
a  =  0 

The  only  restriction  on  h  is  that  it  be  within  the  displacement  assumption 


-Sh  <  h  <  Sh 
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Thus  this  combination  has  an  infinite  number  of  equilibrium  points  and  is 
indeterminate. 

o  Combination  6; 

Sa  <  a  <  ®  and  -Sh  <  h  < 

The  system  equations  are 

0  =  -AM2a 

Ka  a  -  KaSa  =  AM2ra 
Solving  for  a  yields  two  solutions 
a  =  0 

a  =  KaSa/(Ka  -  AM2r) 

This  combination  has  no  solution  and  is  therefore  undefined. 

o  Combination  7; 

-  00  <  a  <  -Sa  and  Sh  <  h  <  00 

The  system  equations  are 

Khh  -  KhSh  =  -AM2a 
Kaa  +  KaSa  =  AM2  ra 

Solving  for  h  and  a  yields 

a  =  -KaSa/(Ka  -  AM2r) 
h  =  -AM2a/Kh  +  Sh 

This  is  a  unique  equilibrium  point  and  the  system  is  determinate. 

o  Combination  8 

-Sa  <  a  <  Sa  and  Sjj  <  h  <  ® 

The  system  equations  are 

Khh  -  K hSh  =  -AM2a 
0  =  AM2ra 

Solving  for  h  and  a  yields 

a  =  0 
h  =  Sh 

This  is  a  unique  equilibrium  point  and  the  system  is  determinate, 
o  Combination  9: 

sa  <a  <  ®  and  Sh  <  h  <  ® 
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The  system  equations  are 


Khh  -  KhSh  =  -AM2« 

Ka  a  -  KaSa  =  AM2ra 
Solving  for  h  and  a  yields 

a  =  KaSa(Ka  -  AM2r) 
h  =  -AM2a/Kh  +  Sh 

This  combination  has  no  equilibrium  point  which  satisfies  the  displacement 
assumptions  and  thus  the  system  is  undefined. 

o  Summary 

For  the  case  of  two  freeplay  springs  with  an  aerodynamic  moment  there 
exist  four  determinate  equilibrium  points,  one  indeterminate  point,  and  four 
undefined  points.  These  results  are  reflected  in  the  data  shown  in  Figure  14, 
Section  3.1.5. 
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BENEFITS  OF  UNIVERSITY/ INDUSTRY  COOPERATIVE  RESEARCH 
ASME  1989  Mechanical  Engineering  Department  Heads  Conference 


Walter  Eversman 
Curators'  Professor 

Mechanical  &  Aerospace  Engineering  and  Engineering  Mechanics  Department 

University  of  Missouri-Rolla 
Rolla,  MO  65401 
(314)  341-4670 


ABSTRACT 


The  University  of  Missouri-Rolla  is  currently  involved  with  McDonnell 
Douglas  Missile  Systems  Company,  St.  Louis,  MO,  in  a  cooperative  research 
program  for  the  Air  Force.  Of  significance  is  the  added  benefit  of 
information  exchange  with  this  type  relationship  --  in  addition  to  the  obvious 
goal  of  doing  good  research  work.  This  university/industry  interaction  has 
developed  during  the  course  of  the  research  activity  and  demonstrates  the 
benefit  of  cooperative  research  programs. 


•  DESCRIPTION 

McDonnell  Douglas  Missile  Systems  Company  is  currently  under  contract  to 
the  Air  Force  Office  of  Scientific  Research  ( AFOSR )  for  a  research  program 
entitled  "Chaotic  Response  of  Aerosurfaces  With  Structural  Nonlinearities." 
This  is  a  joint  program  involving  McDonnell  Douglas  and  the  University  of 
Missouri-Rolla  (UMR).  Dr.  Robert  M.  Laurenson  is  Principal  Investigator  at 
McDonnell  Douglas  and  Dr.  Walter  Eversman  of  the  Mechanical  and  Aerospace 
Enqineerinq  and  Engineering  Mechanics  Department  is  the  Co-Principal 
Investigator  at  UMR.  McDonnell  Douglas  is  doing  nonlinear  aeroelastic 
analysis  development  work  while  UMR  is  designing,  building,  and  testing  wind 
tunnel  models. 

Early  in  the  project,  the  wind  tunnel  model  design  was  used  as  a  project 
for  a  Mechanical  Engineering  senior  design  team  at  UMR.  This  model  was  to 
have  a  rigid  airfoil  with  adjustable  nonlinear  root  support  springs.  The 
students  established  a  configuration  and  fabricated  a  prototype  model. 
McDonnell  Douglas  engineers  attended  the  student's  oral  presentation  of  their 
design  project.  This  allowed  them  to  obtain  a  much  better  understanding  of 
UMR's  design  concept  and  to  provide  suggestions  for  modification  to  the  final 
model  design.  In  addition,  the  students  gained  experience  in  presenting  their 
design  work  to  an  outside  "customer."  The  students'  prototype  has 
subsequently  been  refined  and  has  been  used  in  wind  tunnel  tests  to 
investigate  limit  cycles  which  occur  in  nonlinear  aeroelasticity . 

As  the  result  of  discussions  on  McDonnell  Douglas  analytical  work,  UMR  was 
provided  with  "quick  and  dirty"  flutter  analysis  techniques  McDonnell  Douglas 
uses  during  advanced  design/concept  formulation  studies.  These  methods  have 
proven  useful  in  the  teaching  environment  to  introduce  various  aspects  of  the 
flutter  phenomenon  to  students.  They  have  been  used  as  an  introduction  to  the 
more  rigorous  development  of  dynamic  aeroelastic  theory. 
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